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1. EXECUTIVE SUMMARY 


The objective of this task was to study the origin of the phenomena 
that contribute to the overall noise in photonic phased array systems. In 
particular, the task included the study and the characterization of the 
semiconductor laser low frequency noise and its 1/f behavior. 

The approach adopted for this study included three parts. In the 
first, a comprehensive study of the effect of phase noise from lasers and 
other sources on photonic phased array systems was considered. This study 
considered the influence of additive, as well as multiplicative noise in 
photonic systems, and resulted in the important conclusion that the effect 
of uncorrelated multiplicative phase noise in a phased array system is 
diminished as the number of array elements is increased. The study 
further revealed that, by contrast, the additive phase noise of the array does 
not diminish with the number of elements. Thus for an array with a large 
number of elements, the overall signal-to-noise ratio will be independent of 
the size of the array. Results of this segments of the study pointed to the 
significance of including both the additive and multiplicative noise sources 
in system designs of photonic phased array antennas. 

In the second segment, the influence of modulation on the noise of 
semiconductor lasers was experimentally determined. In particular, the 
dependence of the low-frequency noiseon injection-current modulation in 
external cavity semiconductor lasers was examined. This noise can limit the 
performance of photonic sensors and phased array antenna systems, and 
can reduce the sensitivity of photonic and fiberoptic sensors. 
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The approach followed for the study of the 1/f noise in the spectrum 
of semiconductor lasers was the consequence of the need to examine the 
noise for a single mode laser, based on the findings of the second part of 
the study. This approach included a model of the laser as a noise-driven 
resonant amplifier. A dynamic partial wave model devised for this system 
was then used to obtain the spectrum of the output field of the laser. This 
model produced interesting and suggestive results, including a power- 
independent term in the Schawlow-Townes formula for the linewidth of the 
semiconductor laser, resulting from the multiplicative noise due to electron 
density fluctuations in the gain medium. The model, however, failed to 
directly account for the 1/f behavior in the spectrum as a consequence of 
the laser action. 


2. INTRODUCTION 


Noise is typically the limiting factor in the performance of photonic 
RF systems. The influence of noise in the stability of fiber optic 
distribution links, previously studied in detail at JPL, pointed to the need 
for a comp-rehensive examination of the influence of the phase noise of 
semiconductor lasers in photonic multi-element arrays. The 
Semiconductor Laser Low Frequency Noise Characterization task was 
aimed at addressing this problem. 

The approach used in this task was devised to systematically consider 
the influence of low frequency noise of semiconductor lasers, both 
experimentally and analytically, and then develop a model for the origin of 
the 1/f noise in the laser spectrum. In this final report results of the study 
will be presented. The presentation is divided in three segments. In the 
first part the effect of the phase noise on the performance of phased arrays 
is analyzed. In the second part, the influence of modulation on the laser 
noise is reviewed, and experimental results are presented. The third 
segment includes the presentation of a model for the noise of the 
semiconductor lasers which yields the influence of the multiplicative noise 
on the linewidth manifested as an additional term ‘o the Schawlow-Towncs 
formula. 

While the study did not reveal the origin of 1/f noise in the spectrum 
of semiconductor lasers, it laid important groundwork for future 
examination of this topic, which is difficult to analytically track, yet is of 
fundamental importance in the performance of photonic based RF systems. 
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3 Phase Noise in Photonic rf Arrays 


The function of rf photonic phased arrays is fundamentally based on the 
phase control of individual elements. Thus any noise due to individual com- 
ponents in the system limit the system performance. Basically the total 
noise in. the system is determined by the combination of both additive and 
multiplicative components. The additive term is produced by a variety of 
sources including the thermal noise of dissipative elements, the shot noise of 
the photodetectors, and the laser relative intensity noise (RIN). The additive 
term is readily measured and. in most instances, is easy to calculate; so its 
influence on the performance of the system is generally straightforward. The 
additive noise is always present in a system, and can be determined indepen- 
dent of the presence of any signal. The multiplicative noise originates from 
low frequency gain and path length instabilities, and is only present when 
the signal propagates in the system. This term usually has a 1/P behavior 
with frequency, and thus is typically dominant in frequency regimes close to 
the carrier. 

Our model for the additive and multiplicative noise terms deals with the 
uncorrelatcd componen s. The model is given in detail in Appendix I and 
essentially generalizes the case of noise for a single element with unity gain 
to the cane of a phased array with M elements. The analysis leads to the 
conclusion that for uncorrelated white Gaussian noise, the SNR for an M 
clement array is given by. 


SAT? 
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where SNR is the signabto-noise ratio of the M element array, Eo is the 
amplitude of the input signal, 6<p(t) is the multiplicative phase noise, and 
N(t) is the additive noise term. In the above expression () represent time 
averaging. 

From the cibove equation it is readily concluded that the additive term 
m t he SNR remains constant, independent of the number of elements in the 
array. The multiplicative term, by contrast, does depend on the number of 
arrays, and is reduced by factor M. As M increases, the multiplicative noise 
term reduces to a negligible factor, and the noise performance of the array 
is determined by the additive noise contribution, only. 

As mentioned above, laser RIN is a multiplicative noise source. Thus 
in designing photonic rf arrays the use of the above equation permits the 
determination of the type of the laser best suited for the particular system. 
Since the RIN for the injection current modulated lasers is generally much 
larger than that of the solid state YAG lasers, the combination of the additive 


noise level and the number of elements required will dictate which type will 
produce the highest performance in a given array. 
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4 Influence of Modulation on Semiconduc- 


tor Laser Noise 

In this task, the dependence of the low frequency intensity noise on injection- 
current modulation was studied. The results of an extensive literature search 
were employed to determine the most advantageous manner in which to pro 
eeed. It has been shown previously that low frequency intensity noise (0 * 

1 GHz) in injection modulated semiconductor lasers is upconverted to the 
vicinity of an RF modulation signal. The heart of the problem, therefore, lies 
in determining the origin of the low frequency fluctuations in the laser diode, 
and how the low frequency noise is affected by injection current, modulation. 

Two competing views exist regarding the origin of the low frequency 
fluctuations in semiconductor lasers. One camp contends that low frequency 
intensity fluctuations are due to trapping of carriers in the semiconductor 
medium with a 1/f power law frequency spectrum, while the other view holds 
that competition between longitudinal modes of the laser diode causes an 
enhancement of low frequenc actuations. The initial approach we adopted 
wa*‘? to control the amount of coupling between longitudinal modes of the 
laser, and to observe the effect on the low-frequency fluctuations and close- 
to-carrier RF phase noise. 

To investigate* the effects of mode coupling, an external cavity semicon- 
ductor laser diode was constructed. The external cavity ha* the desirable 
effect of reducing the longitudinal mode frequency separation from approxi- 
mately 450 GHz to 1 GHz, within the range of injection current modulation 


frequencies. Then, the mode-coupling can be affected by injection current 
modulation at appropriate frequencies corresponding to the round trip light 
time of the cavity and its harmonics. 

Testing of the external cavity semiconductor laser was performed under 
various conditions of bias current and modulation. Preliminary results em- 
ploying mode-locking with a strong sinusoidal modulation at 1 GHz caused 
a reduction in the low frequency intensity noise of the laser diode in ex- 
ternal cavity by approximately 10 dB. The cause of the noise reduction is 
presumed to bo the reduction of competition between the lasing modes of 
the laser when mode-locked. However, an unexpected result was obtained 
when the 1/2 sub-harmonic of the mode-locking frequency (i e. , 500 MHz) 
was employed: the low frequency intensity noise was reduced by another 10 
dB for frequencies above 5 kHz. These results were presented at the Third 
Annual DARPA Symposium on Photonics Systems for Antenna Applications 
(PSAA-III), on January 20th 1994 in Monterey California, in a paper enti- 
tled: 'Influence of Modulation on Noise in Semiconductor Lasers," by R.T. 
Logan Jr. and L. Maleki 

5 Multiplicative Noise and Laser Linewidth 

A theoretical model for the mode-competition noise reduction in a multi- 
mode laser was developed ar,d computer simulations were performed. The 
model is based on a Fabry- Perot. cavity containing a medium with a time- 
varying index of refraction. In the simulation, the center frequency and 
amplitude of the dominant mode whs tracked to estimate the frequency and 
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intensity fluctuation power spectra. The computer random-number genera- 
tor was used to simulate index of refraction variations with a white power- 
speetral-density; the model appeared to generate 1/f® -type noise similar to 
that observed in the laser, and proport ional to the noise amplitude. 

This initial theoretical result served as tentative confirmation of the mode- 
competition hypothesis. Theoretically, then, the i/'f noise should disappear 
completely if mode- compet ition was completely suppressed by perfect mode- 
locking. It should also disappear in the rase of a purely single- mode laser. 
However, further searching of the literature, and experience revealed that a 
residual level, of 1/f noise is often observed even in single-mode distributed- 
feedback laser diodes. The computer model was applied to the single-mode 
laser diode a.gain, incorporating index of refraction variations with a w'hite 
power-spectral- density. Again, 1/f- type frequency and intensity fluctuations 
were observed in the simulated output. 

To verify the computer result, an analytic solution was derived for the 
problem of the fluctuations of the transmission of a resonant cavity due to 
variations in the feedback parameters. The analytic model is an extension 
of a hneai* feedback system model that includes the effects of time-variation 
of the feedback components. In this linear model, the transmitted wave at 
any time instant is comprised of a summation of ’partial waves" that have 
existed in the cavity for various lengths of time. In the ca.se of no feedback 
fluctuations, the partial waves may be summed analytically, and the familiar 
Fabry- Perot cavity result is obtained. As the Q of the cavity is increased, 
more partial wave's are stored in the cavity However, when the feedback 
parameter is fluctuating. I he output wave is a summation of many partial 
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waves that have sampled the feedback fluctuations over increasing lengths of 
time. It is seen through this analysis that in a high-Q cavity, the fluctuations 
of the cavity are amplified by the fact that many partial waves make up the 
output wave at any instant, and each partial wave carries a history of the 
cavity fluctuations with it. 

The analytic model results serve to confirm and explain the computer sim- 
ulation results. The fluctuations of the output wave clue to refractive index 
variations with white power spectral density acquire a 1/f 2 power spectral 
density at frequencies above the -3 dB point of the cavity, but remain white 
at frequencies within the cavity bandwidth. This noise enhancement effect 
may explain why high-Q cavities typically have relatively poor long- terra 
stability. 

The linear model is not sufficient to model an oscillator, such as a laser. 
However, it provides interesting new insights into the effects of noise in res- 
onators and resonant amplifiers. In the next phase of work, the time-varying 
feedback and partial wave model was incorporated into a non-linear analysis 
of the resonant cavity above the oscillation threshold. This formalism en- 
tailed modeling the laser as a noise driven resonant, optical amplifier that has 
random gain and phase fluctuations. The model yields an additionaJ power 
independent term in the Schawh w- Townes expression for the linewidth of 
the laser. As the power is increased, the .nodel predicts a re-broadening of 
the linewidth, as generally observed in the case of semiconductor lasers. 

The study also included an analysis of electron density fluctuations in 
the laser gain medium, from which the spectral power density of noise due 
to these fluctuations was obtained. Finally, the results v/erc used to obtain 
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an estimate of the minimum linewidth of the semiconductor laser. 

This segment of the work is the basis for the dissertation of R. T. Logan, 
which will be submitted in the future in partial fulfillment for the Ph. D. de- 
gree in Electrical Engineering, Electro- Physics at the University of Southern 
California. The portion of the dissertation supported in the task is presented 
in Appendix II. 
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APPENDIX I 


Proceedings of SPIt 0/E Lase '94, Optoelectronic Signal Processing for Phased- 
Array Antennas IV, 26-27 January 1994, Los Angeles, CA. 

EFFECTS OF PHASE NOISE 
FROM LASERS AND OTHER SOURCES 
ON PHOTONIC RF PHASED- ARRAYS 

Ronald T. Logan Jr. and Lute Maleki 

Jet Propulsion Laboratory, California Institute of Technology 
4800 Oak Grove Drive, Pasadena, California 91109 

ABSTRACT 

The beam pattern of a linear phased-array antenna system employing a photonic 
feed network is analyzed using a model for the individual feed element noise 
including both additive and multiplicative equivalent noise generators. It is 
shown that uncorrelated multiplicative noise power of the individual feeds is 
reduced by a factor of N in the output of an N-element linear array. However, the 
uncorrelated additive noise of the individual feed paths is not mitigated, and 
therefore will determine the minimum noL. floor of a large phased-array 
antenna. 

1. INTRODUCTION 

In phased-array antennas, the beam pattern depends critically on the phase 
control of the signals at the individual antenna elements. The ability to feed and 
adjust the phase of the microwave signals to the individual radiating elements 
using optical fiber and photonic components offers -ebvious advantages in size, 
weight, mechanical flexibility, and cross-talk, compared to metallic waveguides 
and phase-shifters Various phased-array antenna system architectures with 
photonic feed networks have been proposed, however, the issues of phase 
stability and signal purity are not typically addressed in these proposals. 
Therefore, determination of the acceptable phase noise contribution of the 
individual active feed components has been problematic. In this paper, the 
general factors contributing to the phase stability of an array feed network are 
outlined, with particular attention paid to the type of noise encountered in 
photonic feed elements. It is shown that the analysis of array phase stability must 
consider both additive and multiplicative noise generation processes, and that 
the additive noise of the active feed components will limit the phase stability of a 
large phased-array. 

2. PHASED-ARRAY SYSTEM MODEL 

The general architecture of a phased-array antenna system comprised of M 
elements is depicted schematically in Figure 1. In this analysis, the phase noise 
contribution due to the array feed and antenna elements only is calculated. The 
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effects of the source phase noise will be common to all elements, and may be 
treated in the usual manner for a single antenna element. 



Main Lobe 
Combined 
Output 


N n : Additive Noise Generators 
&fr r : Multiplicative Noise Generators 


Figure 1. Noise Sources in M-Element Phased-Array Feed System 


The feed system is driven by a common source oscillator whose power output is 
divided M ways. The phase delays required to point the antenna beam are 
generated in the separate feeds. To study the effect of noise on the elements, we 
assume a simple case in this model: the signals acquire equal delays in the feeds, 
and are then recombined in a second M-way power combiner. Thus, the output 
signal amplitude is scaled to be equal to the input signal amplitude, to facilitate 
comparison between a single-element antenna system and an array. 

3. NOISE PROPERTIES OF PHOTONIC FEEDS 


The photonic feed elements contain active components such as laser diodes, 
photodiodes, and amplifiers. The phase noise contribution of a microwave fiber 
optic feted system is therefore comprised p* an additive noise term and a 
multiplicative noise term. Laser relative intensity noise (RIN), shot noise, and 
thermal (Johnson) noise are additive noise sources which are present at all times 
independent of signal level. Low frequency gain or path-length instabilities that 
modulate the microwave signal amplitude and phase are multiplicative noise 
sources that are only observed when a signal is present. As shown in Figure 2, 
additive noise usually determines the noise floor at higher offsets from the carrier 

frequency, but multiplicative noise often has a l/f a power law characteristic, 
l<a<2, so is typically dominant close to the microwave carrier frequency [1]. 





Figure 2. Typial Phaie Noise of Photonic Feed Sy*ern at 10 GHz 


Additive noise sources due to thermal and shot processes or laser RIN are 
independent random processes and therefore car. be assumed uncorrelated 
between the feed elements. Multiplicative noise may or may not be uncorrelated 
between elements, depending on its origin. For example, thermal expansion or 
vibration of all the optical fibers in the feed network may produce a phase 
modulation that is common to all elements, whereas, laser or amplifier-induced 
I/f gain fluctuations will be uncorrelated betw een elements. Noise that is 
common to all the elements may be referred to the source srillator and treated 
as if the array were a cingle element. In the analysis that follows, all of the noise 
sources in the individual feed elements are assumed to arise from independent 
random processes, so can be treated as uncorrelated sources. 

It is noted that the multiplicative noise is not detectable by a standard noise 
figure measurement in the microwave signal frequency band. In fact, it is 
difficult in practice to predict the amount of multiplicative noise in an amplifier 
or laser diode, because the noise le_yel may itself be a function of the modulation 
signal frequency or amplitude. Therefore, the amount of multiplicative phase 
noise is usually determined empirically. 

4. ANALYSIS 

We now proceed to calculate the effect of the noise added by the feed elements 
on the total array performance. Consider first a unity gain single feed element 
consisting of a fiber optic link and electronic amplifiers, with a microwave input 

signal E m (0® of constant amplitude E, at microwave frequency 0 ) e . The 
output signal ampli h 'de for a single feed element may be written 

+ N\t) (1) 

where N'(t) represents an additive Gaussian noise process with a white power 
spectral density from dc to well above eo and <5<p(f) is a small multiplicative 
phase noise term. The multiplicative noise term can, in general, be complex and 
thus also represent gain fluctuations. Factoring out the sinusoidal variation at r o, 
yields the slow time-variation of the output field around the microwave carrier 
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+ NV)e"™ 

= E 0 e> m '+N(t) 


( 2 ) 


Limiting the analysis to a band of frequencies So) in width around co e , the 
additive noise term can be written as a random phasor: 
N(f) = N'(f =8r(t)+j8i(t), where |5i(0) = jfirff)), and are assumed to be 
independent Gaussian noise processes with white power spectral density from 
dc to 8(0 (2. Note that as the input signal amplitude E e is decreased to zero in 
Equation (2), the multiplicative noise term vanishes, but the additive noise term 
Nit) is unchanged 

The total array noise is now calculated by using the M individual feed element 
expressions from Equation (2) in the standard calculation [2] of the array output 
field. For a linear array at a steering angle 0, the field distribution in the far field 
of the array as a function of observation angle Q may be expressed as 

E„,,(0,0,t) = + N.<o}. (3) 

Figure 3 is the calculated radiation pattern for a linear array of ten antennas vs. 
observation angle 6 when the steering angle 0 = 0. 



Figure 3. Calculated radiation pattern vs. observation angle for a 10-element linear array. 

We wish to investigate the magnitude of the amplitude and phase fluctuations of 
the main lobe versus the? number of elements in the array. If the noise sources in 
the feeds are urtcorr elated, the statistics of the time-variation of the amplitude 
and phase of the main-lobe peak will be independent of the steering angle <t>. 
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Therefore, for the purposes of this noise analysis, the array can be modeled 
simply as M eqi al-length feed systems sandwiched between two back-to-back 
ideal M-way power splitters, as depicted in Figure 1. Now, the output field 
amplitude is equal to the input field amplitude, independent of the number of 
elements M. 

At the main lobe peak, xcosO - <t>, and the signal amplitude at the output of the 
M-way coupler in Figure 1 corresponds to the main lobe peak amplitude of the 

antenna pattern, divided by VM . For small-angle phase noise <$b(f) « 1 radian, 
the time variation of the output field can be written 



First consider the output of an "array" comprised of only one feed element with 
unity gain. In the above equation, this corresponds to the case of M=l. The 
output field is then given by 

£«*(*) = E„ + jS<t>(t)E c + N(t). (5) 

which is just the sauic as the expression for a single feed given earlier in Equation 
(2), as required. The signal and noise power are proportional to the time-average 
of the squared-magnitude of the output voltage: 

p„, ~ +{jS<pEX ) + (/5*E. (/5b)- E/) + (jS$E c N ' ) (6) 

v +{NE/> + (N(;WE/) + (NN-) 

-where the angle-brackets { ) denote the time-average of the enclosed quantity, 
and the explicit time-dependence of the random functions has been dropped for 
clarity. For independent zero-mean noise processes, the time-average of products 
of the constant and randor ' terms are zero, because we have assumed zero-mean 
random noise processes. Also, the average of the product of two uncorrelated 

terms, such as (N(;50)‘), is zero. But the time-average of the square of any signal 

or noise term is non-zero. Thus, the ratio of signal to noise power (SNR) for one 
feed is 



Similarly, for the case of M parallel feeds with independent equal-amplitude 
multiplicative and additive noise sources, the output field is given by 
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( 8 ) 


The individual powers of equal-amplitude, independent (therefore uncon elated) 
noise sources may be added linearly. Thus, the equivalent noise voltage due to 

the sum of M uncorrelated equal-amplitude noise sources is just Vm times the 
amplitude of a single noise source. Since ail cross-terms between uncorrelated 
noise sources average to zero, we can write the output field for an M-element 
array in terms of a single multiplicative noise source S^(f) and a single additive 
noise source N(t) 



Now, the ratio of signal to noise power in the combined output of an M-element 
array is 



The multiplicative noise power is thus mitigated by a factor of M in the 
combined output of an M-element array, whereas the signal power and additive 
noise power are unchanged from the single-element case. Figure 4 illustrates 
these results for an array of ten elements compared to a single element. 
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Figure 4. Comparison of Singlo-elemcnt and 1 0-elemc nt Phase Noise 

The dose-to-carrier l/f® multiplicative noise of the combined array output is 
reduced by 10 dB compared to the noise of an individual element. The white 
additive noise power further from the carrier is unchanged from the single- 
element case. TT\is behavior is analogous to the improved frequency stability 
obtained from an ensemble of oscillators, compared to the stability of a single 
oscillator. This property may therefore make phased-array antennas more 
desirable than single-element antennas for applications in which high levels of 
long-term phase stability are required. Alternatively, this property relaxes the 
requirements on multiplicative phase noise for the elements of a large phase 
array. It is emphasized that the additive noise requirements are not relaxed, 
however. 
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5. SUMMARY 


It was shown that as the number of array elements M is increased, the effect of 
uncorrelated multiplicative phase noise of the feed elements on the total array 
stability is diminished. However, the uncorrelated additive noise of the feed 
elements is not diminished, so that the signal-to-noise ratio becomes independent 
of array size for large enough M. It is therefore important to quantify both the 
additive and multiplicative noise of the feed elements to correctly predict the 
total array phase stability. 
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Abstract 


A general formalism is developed for analyzing the output field 
fluctuations of a resonant optical amplifier that has random gain and 
phase fluctuations by a modification of the partial wave model. The 
partial wave model provides an intuitive physical picture for the effect 
of multiplicative noise in the optical resonant amplifier. Then, by mod- 
elling a laser as a noise-driven resonant optical amplifier, it is shown 
that multiplicative noise generates an additional power-independent 
term in the Schawlow-Townea formula for the line width of the laser, 
and a re-broadening of the linewidtb at high output power levels, as is 
typically observed in single-mode semiconductor lasers at high power 
levels. Time-varying complex gain constitutes a multiplicative gain 
and phase noise source that is transformed to the output field in a 
different way from additive noise sources, such as spontaneous emis- 
sion. The formalism developed allows calculation of the transmitted 
field power spectra of amplitude and phase fluctuations due to multi- 
plicative noise. It in shown that the effect of muUiplicative noise on 
the transmitted optical field in enhanced as the net round-trip gain in 
the resonator is increased, so that the multiplicative noise ultimately 
determines the minimum Hnewidth for the resonant amplifier in the 
high-gain limit. A detailed analysis of electron density fluctuations in 
a semiconductor laser gain medium is also performed, from which the 
power epectral density of multiplicative phase noise due to electron 
number -density fluctuations in the gain medium is obtained. This re- 
sult is then used in the multiplicative laser noise model to obtain an 
estimate of the minimum linewidth of a aemiconductor laser. Although 


the multiplicative noise analysis is applied to a resonant optical ampli- 
fier, the formalism is general and should be applicable to the analysis 
of other types of feedback systems perturbed by multiplicative noise. 
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1 Introduction and Overview 


A formalism is developed for analyzing the output field fluctuations of a 
resonant optical amplifier with a time-varying complex gain medium, by a 
modified partial-wave model. A linear model of a laser as a noise-driven res- 
onant optical amplifier is then modified using the results of the time-varying 
partial wave analysis. This new laser model provides an intuitive pictUTeof a 
mechanism for power-independent linewidth and linewidth rebroadening in 
semiconductor lasers. 

In a resonant amplifier, time-varying phase and gain constitute multi- 
plicative noise sources whose effects are transformed to the output field in a 
different way from additive noise sources, which are typically treated as ad- 
ditional inputs to the system. The inclusion of multiplicative gain and phase 
fluctuations in a. partial-wave analysis of an optical resonant amplifier is the 
principal novelty of this work. The partial wave model provides a simple 
physical picture for the effect of multiplicative noise in the optical resonant, 
amplifier. The resonant simplifier is modeled as a delay-type feedback system 
(Figure 1) with input , output., and net round-trip loop gain less than unity, 
as shown in Figure 1. The output field at any time instant, is comprised as 
a sum of partial wa ves; the “oldest ' partial waves have travelled more times 
around the loop, so that t he; "memory" of the cavity extends further back in 
time. When no multiplicative noise is present, the width of the cavity reso- 
nance peak decreases as the net gain, and hence the cavity memory time, is 
increased. However, in the presence of multiplicative noise, the older partial 
waves accumulate progressively larger random phase fluctuations, leading to 
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increasing output field fluctuations. For high enough net round-t.rip gain, 
the multiplicative noise-induced fluctuations will dominate the. output field 
fluctuations, and the linewidth will not derease with further increases in net- 
gain. Eventually, as the gain is increased further, the later partial waves 
accumulate so much phase that they begin to destructively interfere with 
the earlier partial waves, leading to large output fluctuations, and increasing 
linewidth. 

By viewing the laser as a noise-driven re^ >nant amplifier, and including 
the results of the partial- wave analysis, a model is obtained for the behavior 
of the laser with multiplicative noise as the net round-trip gain is varied. It 
is shown that multiplicative noise learis t o an additional power- independent 
regime in the Schawlow-Townes formula for the linewidth of a laser at in- 
termediate power levels, followed by linewidth rebroadening.as the power is 
increased further. Firorn this result, an estimated value of 120 Hz for the 
minimum powe-.r-mdopendent linewidth of a typical single-mode DFB semi- 
conductor laser due to electron number-density fluctuations is obtained. 

Multiplicative noise arises in an optical feedback system when a multi- 
plicative factor affecting the field, such as gain or round-trip phase, fluctuates 
in time. Additive noise arises from photons which are added randomly to 
the field, such as by spontaneous emission. In addition to the fundamen- 
tal addtive noine due to spontaneous emission, sources of multiplicative gain 
and phase fluctuation will always be present in the components of a laser 
system at. some level. It. has been shown (12] that statistical fluctuations in 
temperature and density can have significant, effect on the optical phase of 
a wave passing through an optical fiber at. a constant temperature TV For 
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example, in a fiber laser, this source of fluctuation constitutes a fundamental 
source of multiplicative phase noise that arises from the same thermodvamie 
considerations responsible for the unavoidable Johnson noise in electronic- 
components. and is a basic consequence of the Fluctuation-Dissipation The- 
orem [16]. In other laser systems, the net gain and optical length of the laser 
cavity can fluctuate due to dye-stream width fluctuations, mirror vibrations, 
pump power fluctuations, or other technical factors. In typical laser systems, 
it is these technical sources of multiplicative noise that determine the actual 
observed linewidth [1], and not fundamental additive quantum noise due to 
spontaneous emission, as predicted by the well-known modified Schawlow- 
Townes formula [19], [21]. 

The Schawlow-Townes formula for laser linewidth predicts an inverse de- 
pendence on the output power given as 


where 6v c is the ‘'cold cavity” linewidth determined by the losses, h is 
Planck's constant, v is the oscillation frequency, P is the output power, 
and AV A'i are the populations of the upper and low r er states of the atomic 
transition responsible for the optical gain. For large-scale gas lasers with 
high-rcflectivity mirrors, equation (1) predicts linewidths on the order of 
several hertz for output powers in the milliwatt range. In real lasers, the 
observed linewidth typically several orders of magnitude larger than this 
quantum linewidth prediction, due to technical sources of noise. 

A notable exception is the semiconductor laser. The combined effects 


of large output coupling (i.e., low Q cavity) and amplitude-phase coupling 
through the carrier density [2]. [3] yield a predicted fundamental quantum 
linewidth in the kHz to 1GQ MHz range for typical power levels. The pre- 
diction of the appropriately modified linewidth equation (2] for typical facet 
reflectivity of 32 percent yields an estimated linewidth of approximately 30 
MHz at a power level of 10 milliwatts, and experiments are typically in good 
agreement with theory for lower power levels. However, a power-independent 
linewidth of several MHz [11], [22] and linewidth rebroadening (6), [5] in 
single-mode semiconductor lasers is typically observed at high output power, 
which is not. predicted by the modified Schawlow-Townes relation. The excess 
linewidth is detrimental to applications requiring a high degree of coherence, 
such as atomic spectroscopy and coherent communications systems. 

The origin of the excess linewidth in semiconductor lasers is not so clear 
as in the case of large-scale lasers, since the typical macroscopic technical 
sources of noise, such as mirror vibrations, dye stream fluctuations, and the 
like are absent. Early work on the power- independent linewidth attributed 
it to fundamental multiplicative refractive index fluctuations due to electron 
number-density fluctuations [11]. In that work, the transformation of refrac- 
tive index variations fin to fluctuations of the laser frequency Su was derived 
from the phenomenological relation 

iz.ia. m 

u n 

The refractive index fluctuations were included in a root-mean-square sense, 
oerived from thermodynamic consul rations regarding the number-density 
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fluctuations of the electrons contained in the active medium. However, this 
model can not explain the onset of linewidth rebroadening. In later work, 
various other mechanisms have been advanced to explain the origin of the 
power-independent linewidth, such as 1// noise [22]. and spatially-dependent 
temperature and carrier fluctuations [24]. Linewidth rebroadening has been 
explained in terms of excess carrier density in the confinement regions of 
quantum well lasers [5], and spatial hole burning [6]. These or other mech- 
anisms. such as injection-current fluctuations, may be operative in various 
combinations in any particular laser. At present, however, there appears to 
be no uniform explanation of a mechanism responsible for both the power- 
independent linewidth and linewidth rebroadening in single-rnode .semicon- 
ductor lasers. 

The dissertation is organized as follows: In Chapter 2, a discrete-time for- 
mula for the slowly -varying complex envelope of the output field is obtained, 
equation (16). Two regimes of the output fluctuations are examined: small- 
angle and large-angle. The small-angle approximation yields equation (30), 
which has the form of a finite-impulse-response filter acting on the multi- 
plicative and input nose processes to produce the output fluctuat ions. From 
this expression, the: impulse responses and frequency responses of the reso- 
nant amplifier for multiplicative and input fluctuations are obtained, finally 
yielding the output noise power spectrum, equation (46). In the limit of high 
net round-trip gain, this becomes equation (47), and the output fluctuation 
power spectrum in dominated by multiplicative noise. This is the r wehanism 
for the power-independent linewidth. Numerical simulation results for the 
large angle regime illustrate increasing amplitude and phase fluctuations due 
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to the random- walking phase and gain of the partial waves Figures 2 - 
19. With increasing net gain, the probability distribution of the output in- 
tensity in the large angle regime moves to lower values as in 17 - 2\ and the 
power spectrum broadens, (22 - 26 [This is the mechanism for the linewidtb 
rebroadening. 

In Chapter 3, for the small-angle regime, the output field power spectrum 
of the resonant amplifier is found to have a Lorentzian lineshape, that has 
width given by equation (52) for a multiplicative noise source that has white 
power spectrum. This result has the same form as the phenomenological 
result of equation (2). In the large-angle regime, the linewidtb is predicted 
to broaden, which is not predicted by equation (2). The resonant ampli- 
fier analysis is finally applied to the laser by modelling it as a noise-driven 
resonant amplifier. The linewidtb of the laser has a power- independent min- 
imum for the small- angle regime, and rebroadens following the onset of the 
large-angle regime. A minimum linewidth for a typical semiconductor laser 
is estimated to be 120 Hz due to electron density fluctuations. 
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2 General Analysis of Resonant Feedback 
System with Multiplicative Noise 

In this Chapter, the standard linear analysis of a delay-type feedback sys- 
tem is generalized to treat time-varying gain and phase perturbations of the 
components in the feedback system. The gain and phase coefficients multi- 
ply the complex amplitude of the field as it passes through the amplifying 
medium, so that, gain and phase fluctuation; constitute sources of multiplica- 
tive. noise. First, to frame the discussion and define terminology, the basic 
analysis of a linear delay-type feedback system without multiplicative noise 
will be reviewed, and applied to the case of a resonant optical amplifier. The 
gain medium is treated as a linear amplifier operating below saturation. This 
amounts to a review of the standard treatment [4] of a resonant optical ampli- 
fier, from which the familiar Airy formula is obtained. Next, a discrete-time 
formalism, valid for a single oscillation mode at a cavity resonance frequency, 
is developed including multiplicative phase and gain fluctuations of the am- 
plifying medium between the cavity mirrors. 

In the standard analysis of a resonant cavity, the transmitted output 
field is computed from an infinite sum of partially-reflected fields inside the 
cavity, or “partial waves.” Cavity resonances occur at discrete frequencies 
when, all of the part ial waves add constructively For a perfectly stable cavity 
with net round-trip gam less than unity, an analytic form for this summation 
obtains. However, in the presence of multiplicative noise, the partial waves in 
the resonator will have randomly- varying phase and amplitude, so the usual 
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analytic summation is not applicable. Therefore, in the present analysis, the 
output field is computed as a random phasor summation of all the partial 
waves contained in tne cavity at discrete instants of time. The case of an 
ideal noiseless input field is examined first to appreciate the effect of the 
multiplicative noise on the transmitted field. Then the more realistic case of 
an input field with phase and amplitude noise is treated. 

Two regimes are examined: a. small-angle linear regime, and a large- 
angle regime. In the small-angle regime, the fluctuating complex phase of 
the partial waves is less than 0.1 radians about the mean. In this regime, 
expressions for the output field amplitude and phase are analytically derived 
in terms of the power spectra of the multiplicative gam fluctuations and ad- 
ditive input fluctuations, and the net round-trip gain. But it is seen that 
the multiplicative noise is enhanced by the round-trip net gain, and that the 
resonant amplifier system acts like a linear finite-impulse-response filter in 
transforming the multiplicative fluctuations and the input field fluctuations 
to the output field. The impulse responses of the resonant amplifier for both 
multiplicative noise and input field noise are found to decay exponentially, 
leading to Lorentzian terms in the output field power spectrum. The fa- 
miliar filtering characteristics of the Fabry- Perot cavity are obt ained for the 
input additive fluctuations. The response for the multiplicative fluctuations 
is enhanced as the gain is increased, whereas the response for the input field 
fluctuations decreases with increasing gain. For large enough gain, the mul- 
tiplicative noise dominates the output field fluctuations. It is shown that 
the enhancement, of the multiplicative fluctuations limits the output power 
spectrum to a minimum width determined by the net round-trip gain and the 
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multiplicative noise strength, in contrast to the prediction of the standard 
analysis. 

For the case of large-amplitude multiplicative noise exceeding the limits 
of the small-angle approximation, numerical simulations are performed using 
simulated noise time series. In this regime, the power transmission of the 
cavity with multiplicative phase fluctuations grows progressively more erratic 
with increasing net gain. For large enough net-gain, the phase of the output 
field "wraps around” more than 27r radians, causing discontinuous jumps 
of the output phase, and fast, spiky output amplitude fluctuations on time 
scales much faster than the cavity photon lifetime. Also, the probability 
distribution of the output intensity becomes more like that of a thermal 
source, and the power spectrum becomes broadened. In this regime, the 
impulse response is no longer exponentially decaying, but becomes erratic 
and enhanced at large times, resulting in enhancement of the low- frequency 
portion of the output field power spectrum. 

Finally, it is noted that although the primary concern of this analysis is 
the optical resonant amplifier, the result* of this analysis should be applicable 
to noise in other types of feedback systems. 

2.1 Review of Standard Linear Feedback Theory 

To begin the development, the standard analysis of a simple delay-type feed- 
back system with linear, time invariant components will be reviewed, and 
then applied to a resonant optical amplifier similar to the approach found in 
Sicgman [4]. Consider the feedback system model of an optical cavity illus- 


trated in Figure 1. The input, output, and internal states are electric fields 
with sinusoidal time dependence. The input field R incident from t he left of 
Figure 1 is assumed to be a sinusoid of frequency * and amplitude Ar, 

R{v) = A h ( 3 ) 

The output field C(ui) exiting to the right also has the same sinusoidal depen- 
dence. It is desired to compute the transfer characteristics of the feedback 
system from the input state /?(*) to the output stave C(ui). 

The system has an input coupler 1 on the left of Figure 1, and an output 
coupler 2 on the right. These couplers have branching ratio pu for the 
reflected amplitude, and 71.2 for the transmitted amplitude. The couplers 
are assumed to be lossless, so that 2 r p\ ti == 1. The system has complex 
forward gain coefficient (?(>') and reverse gain coefficient H(ui). A crucial 
assumption is that the gains G(w) and H(u) are linear, so that the output 
field can be computed as a vector addition of partial waves. If the gain is 
non-linear, then mixing products will be generated at new frequencies. We 
will assume that the cavity behaves as a. linear resonant amplifier, with the 
output field amplitude linearly related to the input field amplitude. 

Throughout the remainder of the analysis, the explicit sinusoidal time 
dependence of R(ui). C(uj) etc., will be understood, and R, C etc. will be 
used instead. For a pure sinusoidal input field R that has been applied for an 
infinite time, the steady-state output field C is the summation of an infinite 
number of partial waves which may be written 

C = (7i7 a )C?/? + (7,72)(/hrt)C i ^ + (Ti'y2)(^/'a) a G 3 H J « 
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+ (7i 'i2)(PiPi) 4 G 4 H 3 R+-- 

= hr/JCRll + PiP 2 GH + (p l p 2 GH)* + (p l p 2 GH) 3 +-} (4) 

The term outside the square brackets represents the transmission of the input 
and output couplers 7 [ 71 . and the forward gain G for first one-way pass 
through the cavity. The terms inside the brackets represent the successive 
feedback terms which experience a net round-trip gain pip^GH on each trip 
through the cavity. For net gain pip^GH < 1 , the quantity inside the square 
brackets may be analytically summed, so that the output may be written 


C = 


(7l7j )GR 


(5) 


l-intoGHY 

Equation (5) has the form of the Airy formula obtained via the standard 
treatment of a resonant optical amplifier found in optics texts [4]. Identi- 
fying pi,p 2 and 7), 72 as the mirror amplitude reflection and transmission 
coefficients, respectively, and setting G — H = G 0 e.~'™ l,e as the gain and 
phase experienced by a field for a one-way pass through the cavity of length 
l and index of refraction n, with c. the speed of light in vacuum, the familiar 
transmission characteristic for a Fabry- Perot resonant amplifier is obtained: 


t( \ c - 

■’ R 1 -PiPiGfe- 3 "-'*' 


( 6 ) 


The power transmission function is the magnitude-squared of this expression: 


PfM __ (~) 

>.„(«)* " 1 + (.PiPiGi) 1 - 2 HtoCl cc*&jnl/c ) ' ' J 

The resonant Amplifier power transmission function is plotted in Figure 


2 for facet reflectivities P\ =■ m - 0.9 and G 0 = 1. The system is resonant 
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at frequencies corresponding to integer multiples m = 1. 2. • • • of the inverse 
round-trip time through the cavity r - 2nl/c: 


2ttti 




(8) 


Equation (5) is the general form of the input-to-output transfer function 
for a delay-type feedback system as a function of frequency. This expres- 
sion is valid for a system with no multiplicative noise, acting on sinusoidal 
input fields. As the feedback fraction is increased, the widths of the trans- 
mission peaks given by equation (7) decrease. The half-power amplification 
bandwidth of a resonance peak is [4] 


6i; (9) 

p\OlG\ ffT 

So, as the net round-trip gain approaches unity, the bandwidth of the reso- 
nant amplifier approahes zero. 


2.2 Linear Feedback System Analysis with Multi- 
plicative Noise 

In thin section, a general formalism will be developed for deriving the sys- 
tem output field C for the resonant amplifier of Figure 1 when the forward 
and reverse gains are perturbed by multiplicative noise. The multiplicative 
noise may be fundamental or technical in origin, as discussed previously. In 
the general case that C and H ha% r e non- negligible gain and phase fluctua- 
tions. the closed-form summation leading to equation (5) can not be applied. 
Instead, the partial waves must b<? explicitly summed to obtain the output 
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at each desired time instant t, and the amplitude and phase of the output 
field at the optical frequency ut m , will fluctuate in time. This analysis will 
compute the slowly-rvarying complex envelope of the output field at a cavity 
resonance u» m , or C(u> m ,t). 

It is desired to analyze the fluctuations of the output amplitude and phase 
of C{u)m, t) on time scales that are long compared to the period of the optical 
frequency w m . As will be described, this can be conveniently done by moving 
to a discrete- time picture of the system, breaking up time into discrete incre- 
ments equal to the time delay for one round-trip through the cavity, r. Aft< , 
the completion of (his work, a similar approach was discovered in Reference 
[7] in which Z-domaiu techniques from digital systems theory were used to 
obtain the transmission and reflection characteristics of Fabry-Perot etalons 
versus the optical input frequency. However, in that work, no time-variation 
of the feedback coefficients waa considered. Instead, the Z-doraain analysis 
was used to simplify the difficult problem of computing the transmission and 
reflection characteristics of complicated multi-layer structures. In the present 
work, a similar discrete-time approach is used on a simple two-mirror Fabry- 
Perot structure, but extended to the case of time- varying gain coefficients to 
compute the noise properties of the output field. 

For the special case of a feedback system operated near resonance, such 
as a Fabry-Perot optical cavity, the throughput is only appreciable in narrow 
bands of frequencies near the cavity resonances which occur at frequencies 
u/ m , as seen in the previous section. All other frequencies are heavily attenu- 
ated . In many practical applications, such as the 

laser, the fluctuations of the cavity transfer function near a resonance fro- 
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quency are of primary interest, since they determine the long-term stability 
of the laser frequency. Also, the system operates at a single cavity resonance, 
so a full frequency domain analysis is not really necessary, arid i t is sufficient 
to analyze the effect, of the gain coefficient fluctuations on the behavior of 
the cavity near the r esonance frequency of interest. 

Now, the resonant amplifier of Figure 1 will be analyzed near a single 
resonance frequency when the gain coefficients G and H are perturbed by 
noise. To keep the analysis simple and the results transparent, initially it 
will be assumed that the input state R is a pure sinusoid coincident with a 
resonance frequency, u; m , with constant amplitude Ar, as defined in equation 
( 3 ). Later, the more realistic case of an input field with addit ive noise will be 
treated. The input and. output couplers 1 and 2 are assumed to be identical, 
so that pi — Pi = P and 71 = ^ = y. It is assumed that the cavity 
length l is time-invariant, as is the nominal refractive index n. All phase 
or gain fluctuations are described by complex zero-mean random variables 
6g(t) and 6h(t). The time-varying forward and reverse gain coefficients due 
to multiplicative noise are then modeled as 

G{u m ,t) == G 0 e-’ w/e e- , * sKt) = G a G{t). ( 10 ) 

H{u m .t) » H„ = Hn € -inu„t/c (ii) 

where G 0 , ff 0 represent the static forward and reverse loss or gain, 
is the static one-way phase shift added to a signal at frequency passing 
through the cavity. At a loop resonance w m . the static phase shift will satisfy 
the condition nu) M l/c -- 2 tm. radians, with m an integer. bg(t) and 6h{t) are 
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complex random variables representing the time-varying phase and gain at 
frequency u, n . These may be written 

W -W) + %"«) (12) 

6h(t.) = 6h'{t) + iSh"(t) (13) 

where the real parts Sg'(t) and Sh'(t.) represent phase fluctuations, arid the 
imaginary parts and 6h"(t) represent gain fluctuations. The statistical 

description and any correlation of the real and imaginary parts will depend 
in detail on the physics of the gain medium. 

The transformation from a continuous- time representat ion to the discrete- 
time representation is accomplished by considering the static phase shifts 
e -mw m i/e d e i a y elements of duration t/2. as illustrated in Figure 3. and rep- 
resenting the continuous-time gain variations of G(t) and H(t) as a discrete- 
time random sequence of complex values G t and H t with dc gain G n and 
H 0 : 


G t ■■= G 0 e~' Sai = G 0 G t . 

(14) 

$ 

it 

II 

ri, 

> 

1! 

sq 

(15) 


separated by time- intervals r. (Throughout the rest of the analysis, discrete- 
time variables will be denoted with a subscripted time index, e.g., A',, to 
distinguish them from continuous-time quantities: -Y(*)). In this discrete- 
time picture. is assumed to be a pure sinusoid at loop resonance 
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frequency u/ m , with slowly- varying amplitude and phase fluctuations, which 
will be referred to as C ( . 

Now. wc apply this discrete-time definition of the gain to the linear anal- 
ysis of the previous section to obtain an expression for the output field C t . 
At time increment t, the forward gain coefficient is G ( , the reverse gain is H t , 
and the input: state is R t For iiotatioual daritv, it is understood throughout 
that time index t - j refers to time t - jr. since time has been quantized in 
units of the round-trip time in the cavity r. The output state is then written 

+~r*p*GiHt- \Gt-1Ht-2Gt-2Rt-2 + ■ • ■ 

+p i Ht.. l C t . i H,. 2 G l - 2 R < -2+ ■■] ( 16 ) 

ft is assumed that the dc forward and reverse gains and static phase are 
identical, so that G 0 = H 0 . This assumption is equivalent to requiring the 
loss (or gain) and time-of-flight to be equal for both directions of propagation 
in che optical cavity. This is a good assumption for a Fabry- Perot cavity, 
since the light passes through the same physical medium in both directions, 
but may not be valid for other cavity geometries. It is further assumed that 
the gain fluctuation rate is slow compared to the time-of-flight r for a round- 
trip in the cavity. Then, it. is valid to assume the forward and reverse gain 
fluctuations are equal over any round-trip time interval r, so the round-trip 
gam may be written G t H t •= //, 2 . Then, defining the time- varying net round- 
trip gain to be p 2 H} = h\. the output field at time t is given by the infinite 
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senes: 


C t - SHtlR' + K'-xRt- 1 + K t ^K t . 2 Rt-2 + • • •]• (17) 

To analyze the characteristics of the gain coefficient fluctuations on the 
output field, first the case of an ideal noiseless input field will be examined. 
The input state is taken to be a constant R t — R for all times t, so it may 
be fact.orrxi out of equation (17). Then, the output field C t is expressible as 
an infinite sum of partial waves in the cavity at any given time instant t: 


Ct « i 2 H t R[l + #r,_i + (tf,-,tf«_ a ) + (Kt-iKt-iKt-i) + • • •]■ (18) 

This infinite series will be truncated to a finite number of M terms, as 
follows. If the gain fluctuations are small, then the amplitudes of the succes- 
sive terms are determined by increasing powers of the static: net round-trip 
gain coefficient, so that the j - 1st term is proportioned to K ] 0 . When the 
mean value of the net round-trip gain is less than unity, K 0 = p 2 H 2 < 1. the 
infinite sum of K> is equal to 1/(1 — K 0 ), as was used previously to derive the 
Airy formula, if a finite number of terms are used instead, the summation 
of K{ from j = 0 to j — M is given by 


M \ is \1 4 1 

£** = 


(19) 


1 -K q 

The value of M is determined by requiring the difference between the infinite 
sum and the finite surn of A/ + 1 terms to be less than a small fraction t of 
the value of the full infinite sum. This condition is written 
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Then, solving for the required number of terms M yields 

" > £uk- t {n) 

We shall require e < 0 01. or the difference between the finite sum and the 
infinite sum of less than 1 percent. For example, if K 0 = 0 9, then from 
equation (21) we obtain M > 43. The time interval Mr is proportional to 
the memory of the feedback system, and increases as the net round-trip gain 
K approaches unity However, it will be seen later that the value of M just 
derived is only valid in the case of small gain fluctuations. 

The output field at any time t is a random phasor sum of the M + 1 
partial waves. This expression may be cast in a more illuminating form by 
substitut ing the full form of the discrete gain coefficients for the K t terms 
given by equation (15). Doing this and simplifying yields: 

C t = y*H 0 + p^Hl e' 2 ' ( *' , - ,) + p A H* 

p7M f{'* M ( 22 ) 


Each successive term of this expression represents an individual partial wave 
t hat has travelled an increasing number of times through the cavity, and has 
thus experienced the gain and phase fluctuations bh t over an increasing time 
span. For the case of no gain or plia.se fluctuations, (i.e.. Sh ( = 0 for all £,) 
this expression may be summed analytically, and collapses to the well-known 


Airy formula, for the Fabry-Perot transmission derived previously. But in the 
present form. each, partial wave has a different random gain and phase, due 
to the multiplicative complex fluctuations Sh t - 

FYom equation (22) it. is seen that the phases of the successive partial 
waves are cumulative sums of the multiplicative noise process 6h t over in- 
creasing times. If 6h t is a white Gaussian complex noise process with zero 
mean and variance <r^, then the complex phases of the successive terms 
represent a random walk process. The variance of a random-walk process 
grows linearly with time, so the complex phase of the A/th partial wave. 

+ 6h t - 2 H + Sh t . M ), will have variance a\ t = “IXfal. 

There are two regimes which we may identify with respect to the plmse- 
sprc.'iding of the partial waves in the output field summation of equation (22), 
which we shall denote as the small-angle and large-angle regimes. The small- 
angle regime corresponds to relatively weak multiplicative noise and/or low 
net round-trip gain, such that the root-mean-square complex phase of the last 
(i.e., (M 4 l)st) partial wave is less than 0.1. In this case, the exponentials in 
equation (22) may be linearized by the small-angle approximation, and the 
summation for C t may be obtained analytically, as will be shown. 

The large-angle regime corresponds to the case when the the rms com- 
plex phase: variance of the older partial waves exceeds 0.1. In this case, the 
small-angle approximation does not obtain, and the full phasor expression 
of equation (22) must, be used to compute the output field C, at each time 
increment t For a large multiplicative noise variance and/or high net 
round-trip gain h' 0 — 1, the random-walking phase fluctuations of the later 
partial waves may become large enough to cause destructive interference with 
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the earlier partial waves, causing large fluctuations of the output field am- 
plitude and phase. Also, the random-walking gain fluctuations of the later 
partial waves will eventually cause significant, deviations from the exponen- 
tial decay of the partial waves given by K } 0 . This causes the value of M to 
increase, which leads to enhanced low-frequency fluctuations of the output 
field. 

First, the small-angle regime will be analyzed. In this regime, the phase 
fluctuations of the output field due to multiplicative phase noise are en- 
hanced as the net round-trip gain is increased. Then, the large-angle regime 
will be examined by computing the full phasor expression of equation (22) 
for constant multiplicative phase noise variance while progressively increas- 
ing the net. round-trip gain. The transition from the small-angle to large 
angle regimes is marked by the onset of rapid, large fluctuations of the out- 
put field amplitude and phase, and a decrease in the average output power. 
The results obtained will be used to calculate the linewidth of the resonant 
amplifier in each regime. In the small-angle regime, it is found that t he mul- 
tiplicative phase noise leads to an out put- power- independent contribution to 
the linewidth of the output field power spectrum. In the large-angle regime, 
the linewidth increases with increasing net gain. Both of these results are 
counter to the decreasing linewidth with increasing net gain expected from 
equation (9) for the standard analysis without multiplicative noise. 
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2.3 Small-angle multiplicative noise regime: linear 
approx imat ion 

In this section, the case of relatively weak multiplicative noise will be treated 
using a linear approximation. The last term of equation (22) represents the 
'“oldest" partial wave that, has travelled M + 1 round-trips in the cavity. The 
complex phase of this wave is an M + 1-step random walk. The small-angle 
regime is defined as the case when the standard deviation of the complex 
phase for the M + 1st partial wave is less than 0.1, or 


V2Mc h <0.\ ( 23 ) 

where M is defined by equation (21). Then, the small-angle approximation 
for the exponential {e* sc 1 + x) can be used to linearize the expression for 
the fluctuating output field so that equation (22) may be rewritten: 

C t « "i 2 H 0 { 1 +■ i(tfhf))/?[l + pi 2 Hi (1 + 2i(6h t „\)) + p 4 //„ (1 + 2t(<5h,_ 1 + dh t ..i)) 
+ • • • + f?MH™ (1 4- 2 i{6h,.i +Sh t . j + ■ ■ ■ + ^. M ))). (24) 

Substituting K 0 — p 2 Hl Into this expression yields: 


C, - ; 2 H 0 (\4-i(6h t ))R\\ + K u {l + 2i(6h i . ] )) + Kj (1-t- 2>(6/i,_ 1 + ih«-,)) 
+ ••••+■ (1 + 2i(f>h t -i + + • • • + 6f>t- a/ ) )] • (25) 

Expanding this expression will yield many terms of 0{6h J ), like {6h'fbh t .k), 
which are small compared to O(fih), and so may be discarded, leaving 
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C t = -■ 2 H 0 R[(l + Ko + K 2 0 +--+f<") + (l + K 0 + Kl + -- +K?)(i6h t ) 

HK 0 + Kl + K^iSht.i.) + (Kl + K't + • • • + 0(2ifift t . 2 ) 

+ ..-+#rf(2i««-M)]- (26) 

The term proportional to 6h t is due to the phase fluctuation encountered 
in the first ‘‘half round-trip” through the cavity. All other terms represent 
“full round-trips” through the cavity, hence the factors of 2 appear in terms 
5hi_.j. Sh(-2 < ■ ■ ' 

The summations of powers of K a may be performed analytically using 
equation (19). The coefficients in equation (26) are then designated as a*, 
where the coefficient of the first term is a 0 i defined as 

( AJ -f L ; - 1 i __ T/M~\ 

*»-- E ^ = -rir" (27) 

j-0 1 “ A 0 

and subsequent terms a*, k > 1 may be written 


a k = n 0 - E 

j~o 

IzW— 

1 - A'„ 1 - A' u 

1 - K 0 

As ilf — oo. at converges to 


ai< == 



(28) 


(29) 
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Now. substituting r.he coefficients a k into equation (26).. for the present 
special case of a noiseless input R, the output field at any time instant f'is 
seen to be a weighted summation of the multiplicative noise samples h t over 
the time interval t — M to t: 

C* = 't* H 0 R,\q( i + ao'iSht + Q\i26ht~ i + • ■ • + &Mi26ht—\i] 

r M I 

= y 2 H 0 R a fl (l + i6h t ) + .£ a* i26h t - k . (30) 

fc= 1 j 

Finally, the factor 1/(1 — K 0 ) may be factored out of this expression to yield 

C t = (l + i6h t + iKo'PJh^ + iK 2 26h t ~2 + • ■ ■ + iK* 2£/p_„] . 

(31) 

This expression represents the transformation of the mult iplicative phase 
and gain fluctuations to the output field at every time instant t. The first 
term is the static gain of the resonant amplifier as given from the standard 
analysis leading to equation (5). Recall from equation (13) that 6hi is a 
complex-valued random scries whose real and imaginary parts represent the 
time-varying gain and phase, respectively, of the amplifying medium. Phys- 
ically, the; real-valued coefficients a* determine the contribution of the mul- 
tiplicative phase find gain fluctuation at the previous time increment t - k 
to the output at time increment t. It. is apparent from equation (31) that 
for !A'„| < 1. the a* magnitudes decrease exponentially wit.h increasing k. 
Thus, noise events further removed in time from the time increment t have 
a diminishing effect, on the output, state. In this small-angle approximation. 
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noise events occurring at times prior to t - At have no effect on the output, 
which is just an analytic statement of the finite memory time of the system. 

In the analysis thus far. fluctuations of the input field R were not, consid- 
ered. In general, the input R ~ R t may also have a time-varying amplitude 
and phase. To treat this possibility, the input state is now defined in the 
continuous-time picture as a sinusoid with mean frequency coincident with 
a cavity resonance perturbed by amplitude noise 6A(t) and phase noise 

m) 


R(t) = {A R ^6A{t)) e ,(w " t+4 * (t)) . (32) 

Assuming td t < 0.1 radian, the sraaLl-angle approximation can again be 
employed, and the time- varying input field R{t) may be manipulated into 
the discrete-time form: 


R t « (Ah + 6 A t ) (1 - )*’""* 

=• (Ah — iAr60( + 6 At + 16 
~ .4/{ f 1 + i6c>t +• e' a '' n< 

- + /!/?,) 

R, =• R 0 (\ + 6R,). (33) 

where the terra of 0(6A)(Ad) was neglected, since it is small compared to 
the terms of 0{tA) and and the noise sources M and b<p are assumed 

to be uncorrelated 

Now, this expression for the fluctuating input field is inseited into the full 
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form of equation (22) to obtain the output at every time increment t clue to 
both input field fluctuations and multiplicative gain fluctuations: 


C t = y 2 G t Rt + ~t 2 p 2 G t Ht-iG t ..iRt-i 

Wp'GiH^G^H^Gt-iRt-* + • • • 

+p 4 F1 l -. l G t .. l H t .. 2 G t -a^t-2 + • • •] 

= q 2 (7([f?f + Kt-\Rt-\ + Ki-\Ki-iRt-2 + • • •]• (34) 

Using manipulations similar to those used to derive equation (31) (the de- 
tailed algebra is relegated to an Appendix) a simplified expression for the 
complex output field amplitude is obtained: 

C t = + M't + K 0 i26h t - X Kli2fih t ^ + • • • + 

+(1 - K*) (*>Rt + K 0 bR t ~\ + K%6R t -2 + K : o 3 £fl,- 3 + ■ • ■ + 

The powers of K a multiplying the 6R terms may be defined as a set of 
coefficients b *: 


h Kc- (36) 

The full expression for the output field can now be written, more compactly as 
summations ot the two noise processes 6h t and bR t in terms of the coeffic ients 
a ^ and 6^: 
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C t = 7 2 H c R 


a 0 (l •+ bh t ) - 

\k-l / U=0 /. 


This expression repiesents a linear approximation to the full form of equation 
(22) for the slowly- varying complex envelope of the output field C t including 
both input field fluct uations and multiplicative phase and gain fluctuations. 
This linear approximation is valid in the small-angle regime as defined by- 
equation (23). 

Now, the power spectrum of the output field in the small-angle regime 
will be calculated. Equation (37) has the form of a finite-iinpulse-response 
filter operating on “input” processes 6h t and 6R t . For linear systems in 
general, the Fourier t ransform of the time impulse-response yields the output 
frequency response of the system [8]. For an input field noise source modeled 
by a Gaussian-distributed random process, the power spectrum of the output 
field is the product of the input noise power spectrum and the magnitude- 
squared of the system frequency response [?]. So, if the impulse responses, 
net round-trip gain K 0 . and the analytical form of the power spectra of the 
multiplicative fluctuations Si t {f) and the input field fluctuations Sr(/) arc 
known., then the output power spectrum Seif ) may be obtained analytically. 

The impulse-responses of the resonant, amplifier output to a delta-function 
of multiplicative noise 6h t or input noise f>R t are defined to be £CV ( and 6Cr ( , 
respectively The multiplicative impulse-response 6Ch t i& the decaying out- 
put field response given by equation (37) due to a delta-function multiplica- 
tive impulse in 6h t applied to the system at time t — 0: 
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f-c,t = 


y 2 H 0 R 2a t Ut 

2yH 0 R\ 

2^H 0 R 
(l-Ko) 


Kt 


e" rt u t 


(33) 


where the decay constant, and u, is the unit step function, with u t = 0 for 
t < 0 . and u t = 1 for t > 0 , and it is assumed that t takes only integer 
values (i.e., time is normalized to increments of r). The first half-round-trip 
(k = 0 term) was taken to be equal t,o a full round-trip. For M > 100 
corresponding to K 0 > 0.95. this introduces a small error of less than a 1 
percent in the output sum. The response to an input field impulse in 6Rt 
is obtained similarly fis the response of equation (37) to a delta-function 
impulse in 6R t applied to the system at time t - 0 : 


= y 2 H 0 R b, a, 

= ~r 2 H 0 R K[ u t 

= 7 2 H n R c~ vt u t (39) 

The Fourier transforms of the impulse responses 6C\,, and 6Cr it yield the 
corresponding frequency responses F/,(f) and F/r(f), where / is the Fourier 
frequency. Both impulse responses represent decaying exponentials, lending 
to Lorentzian frequency responses given by 

/ 2 \fH 0 R \ 1 

{(\-Kjr) L 
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FhU) = 


(40) 


and 



where — I'/r is the rolloff frequency of the Ixirentzian response. 

The power spectrum of the output field is comprised of three terms: a 
deit<. -function term due to the sinusoidal input R. and terms due to the mul-* 
tiplicativc arid input fluctuations. In the absenc e of multiplicative noise, and 
for a pure sinusoidal input field R at a resonance frequency u > m . the output 
power spectrum Seif) is & pure sinusoid of frequency u m . with magnitude 
given by the constant, term in equation (37). This may be written as 


Sc(f) =- (^H o Ra 0 'fs(f) 

- (P* V</> 


(42) 


where the delta function is defined. a.s <5(0) = 1 and 6(x) = 0, x ^ 0, and 
is centered on frequency u) m . The Fourier frequency / is the offset in V 
from the optical frequency u^. This is equivalent to the power gain i. 
the peak of the resonant amplifier transmission ?. given by equation (7j. 
The noise processes 6h t and fiR t give rise to cuiditional terms in the output 
power spectrum. It is assumed that the noise processes 6h t and 6R t are 
uncorrelated, so the power spectra due to each process may be added. The 
term due to multiplicative noise is the product of the multiplicative noise 
power spectral density S^(f) and the magnitude-squared of the frequency 
response F^(f). Similarly, the term due to the input noise is the product of 
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the input field, noise power spectral density 5 r(/) and the magnitude-squared 
of the frequency response Fn{f). The addition of equation (42) and the two 
noise terms yields the total output power spectrum: 


s c(f)= <H/)j + IW)I J W) + IW)| 2 S*(/) (43) 

Substituting for the frequency responses F h (f) and Fn(f) and subsequent 
simplification yields 



Consider the last term of this expression due to the input noise 6R. For 
low gain and/or low multiplicative noise, the signal-to-noise ratio inside the 
amplifier bandwidth will be the same as the inmd signal- to- noise ratio. Phys- 
ically, the last term implies that the output field phase and amplitude exactly 
follows the input: field phase and amplitude fluctuations for rates less than 
the cavity bandwidth uv, which is intuitively correct. For fluctuation rates 
faster than u/ c . the input field fluctuations are rolled off as l// 2 , as expected, 
since they lie outside the bandpass of the cavity As the net gam K 0 is 
increased, u/ c decreases, and the integrated noise power due to the input fluc- 
tuations decreases, implying a decreased spectral width of the output field. 
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Sn(f) 

(44) 


The resonant amplifier therefore produces a narrowband output signal of 
high spectral purity from a broadband noise input when the net round-trip 
gain K 0 is high. This is the expected behavior of a resonant amplifier as 
obtained from the standard analysis leading to equation (7). and is the basis 
of the linear model of the las-. r as a noise driven resonant amplifier [4], which 
we will return to later. 

Now consider the second term of equation. (44) . It is seen that t he second 
term due to multiplicative noise 6h t is multiplied by the factor 1/(1 - K a ) 2 
compared to the last term due to the input field fluctuations 6R t . What 
is significant about this is that as the net gain K 0 is increased, given fixed 
levels of input fluctuations Suif) and multiplicative fluctuations S^if), even- 
tually the output field power spectrum Seif) will become dominated by the 
multiplicative noise. After this happens, the signal- to- noise ratio for Fourier 
frequencies less than a> c will decrease with further increases in the net gain. 

Making use of the fact that for K„ « 1, P -- \fn(K 0 )\ % (1 - K 0 ), the 
rolloff frequency may be approximated as as (1 - K 0 )fr. Making this 
substitution in the second term of equation (44) yields 




hJJjcR 

\1 - Ko 


+ -—~rlry S M 


1 


(2r 
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In the limit as h' 0 — 1. the second term due to multiplicative noise over- 
whelms the last term, so the output power spectrum becomes 
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w - {r=ff (w) + (5^7? s ‘">) < 46 > 

In the limit as K 0 — * 1 , the integrated multiplicative noise power given by 
the second term in equation (45) approaches a constant value, which implies 
that the spectral width of the output signal also becomes constant. This is in 
marked contrast to the prediction of the standard analysis [4] of equation (7), 
which predicts decreasing linewidth with increasing gain. In summary, for 
low net gain when the input field fluctuations are dominant, the width of the 
output power spectrum decreases with increasing gain. But as the net gain 
K 0 — > 1, the multiplicative fluctuations eventually become dominant, and 
the power spectrum of the resonant simplifier approaches a constant value. 
Therefore, a minimum linewidth should be expected for resonant amplifiers 
operated at high gain in the presence of multiplicative noise. To the extent 
that, the multiplicative fluctuations arise from fundamental thermodynamic 
processes [12], this minimum linewidth will be a fundamental limit. 

2.4 Large-angle multiplicative noise: numerical anal- 
ysis 

In this section, the full phasor form of the output field given by equation 
(22) is analyzed in the large-angle regime, by numerical simulation using 
computer-generated random time seiies to represent the multiplicative noise. 
This approach is taken to investigate the behavior of the output field ampli- 
tude and phase when the variance of the complex phase of the older partial 
waves exceeds 0.1 radians, so the small-angle approximation which permits 
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the .analytic derivation of equation (23) no longer obtains. The output field 
time series is computed for a noiseless input field R of unit amplitude at fixed 
cavity resonance frequency u> m , for various values of the net round-trip gain. 
The calculations are performed for a 30()^m-long Fabry-Perot semiconduc- 
tor !as«?r cavity, at 1.3/rm wavelength, with equal facet power reflectivities 
p 2 = 0.32. 

The real and imaginary parts of 6h t are assumed to be correlated, as is 
the case for phase and gain fluctuations arising from electron-number fluc- 
tuations in a semiconductor laser gain medium [10]. The real and imaginary 
parts of 6h t are written in terms of a common noise source <5 tj ( representing 
fluctuations of the refractive index: 

f.h[t) = 6h'(t) + i6h"{t) 

* (' + i) (¥) Mt) im 

The constant p relates the changes in the real and imaginary parts of the 
refractive index due to electron density fluctuations, and is a measure of 
the amount of amplitude-phase coupling for waves in the gain medium. It 
was shown [2], [3] that P is responsible for broadening of the semiconductor 
laser linewidth above the Schawlow-Townes prediction, and so /? is typically 
referred to as the linewidth enhancement factor The refractive index in a 
semiconductor laser gain medium is typically modeled as a linear function of 
the carrier density [?]. Thus, changes in the electron density due to injection 
current noise, thermal fluctuations, or electron-hole recombination processes 
cause correlated changes in the gain and phase of the amplifying medium via 
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Computer simulations were run for increasing net round-trip gains K 0 « 
0.9, 0.95, 0.99, 0.995, and 0.999, which represents three orders of magnitude 
change in the net gain. Fbr each value of K„ an -element complex vector 
of the fluctuating output field Ct is produced. For these net gain values, the 
total numbers M of partial-waves used in the Ruminations were, respectively, 
Af >■> 50, 100, 600, 1000, and 5000. The complex multiplicative noise process 
6h% is simulated by a computer-generated Gauaslan-diatributsd zero-mean 
random time series with standard deviation of ov » 0.1. This implies that 
the system is In the large angle regime for all of the net gain values tested. 
The computer produces a series of random values from a Gaussian distribu- 
tion such that the rms value of an Infinite number of such values would be 
<Th>- 

The computed results for net round-trip gain of K c • 0.9, 0.95, 0.99, 
0.995, and 0.999 are now discussed. To see the qualitative effect of the 
increasing net gain on the output, tbs real and imaginary parts of the output 
electric field C, ere displayed parametrically In the complex plane in pig 2 

- 6, for increasing values of K t . In these plots, each of the N dote on the 
complex plane repre s ents the computed position of the tip of the output field 
vector C t . Figure 2 illustrates the case of net round-trip gain K e » 0.9. The 
tip of the electric field vector describes an arc with marimmw phase excursion 
cf approximately +/• 0.3 radians, but with relatively constant amplitude. As 
the gain is increased to K t ■« 0.95, 0.99, 0.995, and 0 999, illustrate din Fig 2 

- 6, the maximum phase excursion increases and the amplitude exhibits 
progressively larger fluctuations which appear as a swirling pattsrn in the 
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complex plane. If there were no multiplicative noise, all of the dots would 
fall on one point on the real axis, determined by the static gain 1/(1 - K a ). 
The email-angle regime treated previously represents small excursions about 
this point. 

As the output field vector swirls around the complex plane, Figure 4e, 
the magnitude and phase change drastically. The phase and instantaneous 
intensity (magnitude-squared of the electric field) versus time for each value 
of net gain are plotted in Figures i - 16. As the net gain increases, the 
period and amplitude of both the intensity and phase fluctuations increase. 
At the lower gain values, the mean value of the intensity fluctuations is 
close to the value expected without multiplicative noise. However, as the 
gain is increased, the amplitude of the intensity fluctuations increase greatly. 
The intensity becomes more deeply modulated, and the mean value of the 
intensity moves to lower values. This behavior is more apparent from the 
probability distribut ion of the intensity, discussed below. As the net gain is 
increased, the phase fluctuation amplitude also increases ( Fig 12-1 6 ) until 
eventually, the pliaae “wraps around" more than n r adians 

The probability distribution of the power versus increasing net gain is 
plotted in Fig 17 - 2;iAt K„ = 0.9, the probability distribution has devel- 
oped a “tail" stretching to lower values of intensity. As the net gain increases, 
this tail becomes more pronounced, and the mean value of the probability 
distribution decreases. At high gain of K a - 0.999, the intensity probability 
distribution is peaked near zero, nearly resembling the probability distribu- 
tion of a thermal emission source, or a below-threshold laser. The shift of the 
mean intensity to lower values is in qualitative agreement with a theoretical 
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Fokker-Planck analysis and experimental results obtained for a helium-neon 
laser with multiplicative loss noise [15]. This shift of the probability distribu- 
tion to lower values with increasing gain is contrary to the behavior expected 
for the resonant amplifier without multiplicative noise. In the context of laser 
operation, this result implies that the mean value of the intensity probability 
distribution will initially grow as the gain is first increased, but eventually 
will reach a maximum value, and then begin to decrease. 

Inricj 22 - 26the power spectrum of the output field versus increasing 
net gain is plotted. The noise Level progressively increases as the net gain 
is increased. This is not unexpected, given the increasingly large intensity 
and phase fluctuations of the output field as the net gain is increased. In 
the absence of multiplicative noise, the output power spectrum would be a 
delta-function, since no input noise source is included in the simulations. 

It is noted that in all of t hese results, the variance of the multiplicative 
noise time-series was kept constant, and only the net gain was changed. While 
such a large value of multiplicative noise may not typically encountered in 
practice, these simulation results serve to illustrate the effects of increasing 
net gain on the output field. By using a relatively large value of multiplicative 
noise, it was possible to study the qualitative behavior of the output versus 
net. gain using a manageable number of partial waves. In a semiconductor 
laser, net gains of K 0 » 1 - 10“ 5 are typical, requiring approximately one 
million partial waves per time increment to simulate. To obtain a reasonable 
estimate of the output power spectrum, at least 1024 time increments need 
to be computed, which becomes computationally prohibitive. 
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2.5 Detailed description of simulation algorithm 

For each value of K 0 , the details of the simulation algorithm are as follows: 
First, a real-valued random noise vector 6n t of length 2N is generated. Then 
the complex vector 6h t of length 2.V is generated according to equation (47). 
The cumulative sum of 6h t is then computed, resulting in another vector of 
length 2N , that represents a complex random walk process. The indices of 
this vector axe reversed, so that random-walk begins with the last (2A r th) 
element and moves to earlier times. This vector is called iv. Then, the 
output field Ct is computed for each of the N time increments t by summing 
the M -f 1 partial waves whose complex phases are given by the values of 
w. Thus, the first sample of C t is computed by starting at the middle of 
the vector ui (i.e., wjv-i)- and using the previous M values. The complex 
phase of the first partial wave at time increment t is given by u^+.v-i - 
ti/'e + (,v_ 2 , and the last partial wave has complex j ^ase w t +.v-i - 
Then, the magnitude-squared and phase of C t are computed. The magnitude- 
squared of C t is proportional to the intensity of the output field at time t. 
The probability distribution ( histogram ) of the intensity is also computed. 
Finally, a Hanning (raised-cosine) windowing function is applied to the time 
series Ct and a complex fast-Fourier-transform (FFT) is performed on this 
time series. The magnitude squared of the FFT is the computed power 
spectrum of the output field. 
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responding frequency responses Fh(f) and Fr(/), where / is the Fourier frequency. 
Both impulse responses represent decaying exponentials, leading to Lorentzian fre- 
quency responses. The frequency response of the cavity to the multiplicative noise 
process is 



and the response to the input field fluctuations is 
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(2.41) 


where = V/r is the rolloff frequency of the Lorentzian response, and the Fourier 
frequency f is the offset in Hz from the optical frequency u f)i . 

The power spectrum of t he output field fluctuations is comprised of two terms 
due to the multiplicative and input fluctuations. The carrier at the cavity resonance 
frequency w m is shifted to DC due to the discrete-time analysis, and has amplitude 


Sc( 0) = 




(2.42) 


This is equivalent to the output power at the peak of the resonant amplifier trans- 
mission as derived previously in equation (2.5). 

The noise processes 6h t and SR t give rise to complex- valued terms in the out- 
put field. It is assumed that the noise processes 6h t and 6R t arc uncorrcJated, so 
the power spectra due to each process may be added The term due to multiplica- 
tive noise is the product of the multiplicative noise complex power spectral density 
Sh(f) and the magnitude-squared of the frequency response F K {f). Similarly, the 
contribution due to the input noise is the product of the input field noise complex 
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power spectral density Sr(/) and the magnitude-squared of the frequency response 
//*(/). Therefore, the total output field noise power spectrum due to both input 
field fluctuations and multiplicative gain and phase fluctuations is: 


ScU ) = i W)f 5fi(/) 4 IW)| J SkU) (2.43) 

Substituting equations (2.40) and (2.41) for the frequency responses l\(f) and 
Fr{/)s and using the approximation when K 0 ~ 1 that I’ = |^n(A' a )| 1 K u . 

yields the total output field power spectrum: 
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(2,14) 


This is the general expression for the output, field fluctuation power sport rum of a 
resonant amplifier including input field noise and multiplicative noise. 

Consider now the first term of this expression due to the input noise tR Foi 
low gain and/or low- level multiplicative noise, the signal- to- noise ratio for Fourier 
frequencies less than u> c will be the same as the input signal-to- noise ratio Physically, 
this implies that the output field phase and amplitude exactly follow the input field 
phase and amplitude fluctuations for fluct uation rates loss t han t he cavil \ bandwidth 
u; c . which is intuitively expected. For fluctuation rates faster than the input field 
fluctuations are filtered by the cavity, and so are rolled off as 1 j f 2 . As the net gain 
A' 0 is increased, the cavity corner frequency u; c decreases, and the integrated output 
noise power due to the input fluctuations decreases, implying a decreased spectral 
width of the output field. The resonant amplifier therefore produces a nanowband 
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output signal of increasing spectia) purity from a broadband noise input as the net 


round trip gain K 0 is increased. This is the expected behavior of a resonanl amplifier 
as obtained from the standard analysis leading to equation (2.5), and is the basis of 
the linear model of the laser as a noise-driven resonant amplifier, to which we will 
return later. 

Now consider the second term of equation (2,44) . This term due to multiplicative 
noise is new, and not. obtained from the standard analysis of the resonant amplifier. 
The multiplicative noise term is multiplied by the factor 1/(1 - AV) 2 compared to 
the first term due to the input field fluctuations. It is seen that as the net gain 
K 0 is increased, given fixed levels of input fluctuations S'n(f) and multiplicative 
fluctuations 5u/), eventually the output field power spectrum Seif) will become 
dominated by the multiplicative noise term. After this happens, the signal-to-noise 
ratio for Fourier frequencies less than t j r will decrease with further increases in the 
net gain. 

Again making use of the fact that for Ah fcs 1, I' jfnf A' 0 )| ^ (1 - A' t> ), the roll- 
ofl frequency may be approximated a^ -lv (1 — A\)/ r - Making this substitution 
in the second term of equation (2.44 ) yields 



In the limit as I\ 0 -■* 1. the corner frequency - > 0. and the second term due to 
multiplicative noise overwhelms the first term. The output power sped rum then 
becomes 
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fluctuations and increasing linewidth make the large-angle regime somet hing 
to be avoided, .so that a prediction of the onset of the large angle regime is 
sufficient for most applications. 

By the Wiener-Khinchin Theorem, the power spectrum of a stationary 
noise process is the Fourier transform of its auto-correlation function [9]. For 
the case of small-angle multiplicative phase noise, the amplitude fluctuations 
of the output field are small compared to the phase fluctuations, so the field 
auto-corrchition function is determined primarily by the phase fluctuations. 
In the large- angle regime, large amplitude and phase fluctuations occur, as 
seen in the previous chapter. In this regime, the field autocorrelation func- 
tion and linewidth are determined primarily by the amplitude fluctuations. 

First., the amount of linewidth enhancement is calculated in the small- 
angle regime from the power spectral density of phase fluctuations obtained 
from the small-angle analysis. The small-angle fluctuations are found to lead 
to an oul put-power-independent Lorentzian power spectrum for the resonant, 
amplifier. The power spectrum in the large-angle regime broadens with in- 
creasing net gain as the numerical results of the last chapter illustrated. The 
onset of the large-angle regime is derived from the net round-trip gain and 
the- multiplicative- noise variance. 
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3 Multiplicative noise in laser linewidth 


3.1 Analysis of effect of small-angle multiplica- 
tive noise on semiconductor resonant amplifier 
linewidth 

As discussed in the previous chapter, the multiplicative phase fluctuations 
in a semiconductor laser gain medium arc larger than the gain fluctuations 
by the factor 0. So, in the small-angle regime, the amplitude fluctuations 
of the field are small compared to the phase fluctuations. This is illustrated 
graphically in the complex plane in Figure 4a of the previous chapter. From 
the numerical simulation results at the onset of the large-angle regime, it is 
seen that the phase fluctuations are larger than the amplitude fluctuations 
of the field. This means that the multiplicative noise contribution to the 
output power spectrum $c(f) is primarily phase noise for a semiconductor 
laser amplifier in the email-angle regime. 

The output field power spectrum is the Fourier transform of the field 
autocorrelation function. Following the development of Peterfnann [10], the 
output field power spectrum may be derived in terms r-t the frequency fluctu- 
ations of the field when the amplitude fluctuations are negligible. In this case, 
the field autocorrelation function may be obtained in terms of the frequency 
fluctuations of the output field [10] (See Appendix for detailed derivation): 

(C(t)C’(t - T )) - (P)c ( ‘<* >T) eH< A#,>) (48) 

where P is the average output intensity, (4>) is the mean frequency offset 
from the mode frequency u/ m , and the mean square phase fluctuation ( Ad> ;1 ) 
is related to the power spectTum of the frequency fluctuations 5„(/) as 
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(x/T) 2 

The power density spectrum of the field is then obtained as 


( 49 ) 


W c [u>) - (C(f)C’(t - T))e~' { ~ i d,T. (50) 

J -TO 

These expressions relating the power density spectrum of the field to the 
frequency fluctuation power spectrum are valid as long as the the phase 
fluctuations exhibit a. Gaussian probability distribution. The derivation of 
these results is detailed in the Appendix. 

When the multiplicative noise contribution to the output power spec- 
trum Seif ) (derived in the previous chapter) is predominantly phase noise, 
the power spectrum of instantaneous frequency fluctuations is related to the 
power spectrum, of phase fluctuations as [18] 


ScJ.f) = f 2 S c U). ( 51 ) 

Substituting the previous result from the small-angle analysis for the power 
spectrum of output phase fluctuations, the power spectrum of the output 
frequency fluctuations for net gain K 0 approaching unity may be written 

ScM) - f 2 - { —jzy ShV) (52) 

The factors of / 2 cancel, leaving 

ScM) « >$h(f) ( 53 ) 

(wr r 
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The frequency fluctuation power spectrum thus has the same form as the 
power spectrum of multiplicative fluctuations S/,(/). 

Consider the special case when the power spectrum of multiplicative fluc- 
tuations S h (J) is white. Equation (53) then implies that the power spectrum 
of output frequency fluctuations, ScAf)> is 4180 white. In this case, the 
autocorrelation function of the field amplitude is (See Appendix) 


(C(t)C-(t - T)) = (54) 

where t c — 2/ Sc „(0) . The Fourier transform of this expression yields a 
Lorcnl zian-shapod spectrum: 


W C A) = 


2t-(P) 


(55) 


1 + ([u - ( 6 )} fc) 2 
which is centered around the optical frequency + (ri), and has spectral 
width 


miiif 
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Sh( o). 


(56) 


In the small-angle regime, the total linewidth has the functional form: 


tv\ to , ~ *l> 4- i>l'\rn,jt ( 5 ") 

where the first term on the right is the linewidth obtained from the standard 
analysis without multiplicative noise [4]. 

In the standard analysis of the resonant amplifier without multiplicative 
noise [4; discussed in the previous chapter, the linewidth is give by: 
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For K 0 fa 1. this may be written 


1 - K a 


Some other useful relationships are: 


2r _ T _u c _ ( jj^A'ok 

2 xr xr x 2xnl 


The total linewidth of the resonant amplifer can finally be expressed in 
terms of the net gain K 0 as the sum of the standard and multiplicative 
contributions 


1 — 


(1 - K*) 


7 --T,; S, : {0) 

(2 7T)'V 


when K v m 1. The first term from the standard analysis decreases with 
increasing net gain A'<>, whereas, the second term due to multiplicative noise 
is approaches a constant jis K n — ■ 1. 

A noise source that has a white power spectrum up to irifinfc frequency 
has a delta- function auto- correlation function, which is a useful mathemati- 
cal model, but is not physically realistic A more realistic model for a noise 
process has an exponential tiine-corellation proportional to with cor- 
relation constant <>. The Fourier transform of an cxponcntiallv correlated 
noise source yields a. Lorentzian power spectral density, which is white up to 
a Fourier frequency of approximately / = a. and then rolls off as 1 / f 2 for 
higher frequencies. Then, for maetical purposes, a white noise source may 
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be modeled as having a a larger than all Fourier frequencies of interest. In 
practice, the noise may be considered to be white if a is at least as large 
as the cavity mode spacing equation (56) may be used to estimate the 
Iinewidth. 

When a is small, then the approximation of equation (56) is no longer 
valid, and the Iineshape of the field power spectrum will not necessarily be 
Lorentziao. The Iineshape must then be computed from the Fourier trans- 
form of the auto-correlation function of the electric field. The electric field 
auto-correlation function must be first be derived in terms of the frequency 
fluctuation power spectrum, as is done in reference [10]. 

3.2 Onset of the large-angle regime 

From the numerical analysis of the previous chapter, it is seen that when the 
phase of the later partial waves exceeds the small-angle regime of 0.1 radians, 
the output amplitude and phase fluctuations increase dramatically, (Figures 
2 to 16, previous chapter). In this regime, the iinewidth of the resonant 
amplifier can no longer detemined from the phase fluctuations of the output 
field, as in the small-angle regime. As illustrated in Figures 8a-e, the power 
spectrum of the electric field becomes progressively noisier as the net gain is 
increased, corresponding to increased Iinewidth. 

The large amplitude and phase fluctuation* make this regime of operation 
undesirable, so avoiding the onset of this regime is of great practical interest. 
From the pnsceeding numerical analysis, the onset of the large amplitude 
fluctuations just, begins to occur for net gain of K 0 « 0 9 and when the 
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multiplicative noise standard deviation is C\ = 0.1 radians, as illustrated in 
Figures 4a and 5a. For this value of net gain, the number of partial waves 
to achieve 1 percent convergence of the output partial wave summation is 
computed to t* approximately 45. So, the variance of the phase for the last 
(i.e., Mth) partial wave at the onset of the large angle regime is 

y/2Mo h « v / 90 (0.1) « lradian. (62) 

So, to avoid rebroadening of t he linewidth, the net gain should be kept below 
the value that causes the standard deviation of the A/th partial wave to 
exceed the 0.1 radians, i.e., within the limits of the small-angle model, as 
discussed in the last chapter. 

3.3 Application of multiplicative noise model to laser 
linewidth 

In this section, the results of the small angle analysis for the effect of mul- 
tiplicative noise in a resonant amplifier will be applied to a linear model of 
a laser to derive an analytic expression for the linewidth. The laser may 
be considered to be a resonant amplifier driven by an equivalent input noise 
source corresponding to the spontaneous emission in the active medium [Id], 
& [20]. In this model, the presence of the constant amplitude spontaneous 

emission noise input field causes the net round-trip gain to saturate at a 
value slightly less than unity, resulting in a finite linewidth for the resonant 
amplifier. The ratio of the output intensity to the input noise intensity due 
to spontaneous emission determines the net round-trip gain. The wide-band 
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spontaneous emission noise is amplified within the narrow bandpass of the 
resonant amplifier, yielding the laser output. Thus, for high gain, the laser 
output field may be viewed as a wave of constant amplitude and frequency, on 
which small statistical amplitude and phase fluctuations due to the sponta- 
neously emitted photons are superimposed. The modified Schawlow-Townes 
model [21] takes into account the damping effect of the saturated gain above 
threshold on the amplitude fluctuations, which results in a factor of two de- 
crease in the linewidth above threshold. Although the linear laser model does 
not explicitly include damping of the amplitude fluctuations, its simplicity 
makes it easy to see the basic physics responsible for the linewidth and the 
results for the multiplicative noise are easily obtained. 

In both the original and the modified Schawlow-Townes models, the 
linewidth at a given output pow f er level is seen to arise from the additive 
noise source of the spontaneously emitted photons which add random am- 
plitude and phase components to the output field. The power output of 
the laser increases as the saturation level of the gain medium (i.e., pump- 
ing) is increased, due? to increased stimulated emission. However, since the 
spontaneous emission level is proportional to the gain, which is essentially 
unity above threshold, so it remains clamped at its saturated threshold value. 
Thus, the additive random phase contribution of the spontaneous photons 
becomes less significant os the power increases, resulting in an increasing 
signal-to~noi.se ratio and decreasing linewidth Thus, by the the Schawlow- 
Townes model, it is predicted that the linewidth of a laser can he made 
arbitrarily small by operating the laser at an arbitrarily high optical power 
level. We now' derive the linear laser linewidth model, including the effect of 
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multiplicative phase fluctuations. 

Following the development of Siegman [4j. the output power of the laser 
P is related to the spontaneous emission input noise power P, p as 

1 -RG P ~ 2T~ (63) 

where R is the power reflection coefficient of the cavity mirrors, so that 
R = p 2 , T — o 2 = 1 - R is the power transmission coefficient, and G p is the 
power gain per round trip. The quantity 1 - RG P is defined as the “saturation 
amplification 1 ’ of the resonant amplifier . The spontaneous emission power per 
cavity mode can he written 

(f,4) 

where h is Planck’s constant, v is the frequency of the light, c is the speed 
of light in vacuum, n is the average value of the refractive index, / is the 
length of the cavity, and and jVj are the population of the upper and 
lower atomic levels. For low-loss mirrors, R % l, and above lasing threshold 
the round-trip gain RG V % 1, so (G p - 1) (1 - R). Now, substituting this 

and equation. (64) into equation (63), and assuming lossless mirrors, so that 
T = 1 — R. obtains 


From the development of the previous section, the linewidth of a resonant, 
amplifier with gain is given by equation (b\) 


= l -^ + ~r, 5.(0). (66) 

The round-trip power gain RG P is approximately 1. as is K 0 , so we may 
approximate 1 - RG P » 1 - K 0 . Making these substitutions for the saturated 
gain of the laser into the last expression then yields 


x fed -R)\\ hu 1 

- (“SSrJ 2:1 PAW r + jj^SMO) 

- + ^5 SMO) ( 07 ) 

for the linewidth as a function of the output power, where Si\. is defined as the 
“cold-cavity" linewidth, or the linewidth of the cavity with the gain medium 
removed. The first term of this expression is equivalent to .the original for- 
mula derived by Scliawlow and Townes for the linewidth of the laser, which 
is correct up to threshold, and displays the familiar 1/P dependence. Above 
threshold, this term is reduced by a factor of two due to the damping of the 
amplitude fluctuations by the saturated gain [21]. The second term repre- 
sents the multiplicative phase noise of the laser, and is seen to be independent 
of the output power. This equation predicts t hat the laser linewidth should 
initially decrease linearly with increasing output power, but will finally reach 
a minimum value; proportional to the multiplicative noise strength. If the 
gain is increased yet further, the linewidth begin to ‘rebrnaden” at the onset 
of t he large-angle regime. A minimum linewidth and linewidth rebroadening 
is typically observed in single-frequency laser diodes at high power levels, 
corresponding to high net round-trip gain. Tins model has the same power- 


independent character as the linewidth derived by Welford and Mooradian 
•11] using the phenomenological expression of equation (2). But the formal- 
ism used to derive it also reveals that same mechanism can be responsible 
for both a power-independent iinewidth and linewidth rebroadening, as is 
typically observed in single-mode semiconductor lasers. 

It is well-known that semiconductor lasers exhibit a power-independent 
linewidth deviation from the modified Schawlow-Townes formula [11], [22]. 
As an application of this model, the minimum linewidth of a typical single- 
mode semiconductor laser will be estimated due to the effects of electron den- 
sity fluctuations as derived in the next chapter. Consider a single- mode GaAs 
semiconductor laser amplifier of length l = 250/im, which implies round-trip 
time r as 5 picoseconds. Assume multiplicative phase fluctuations due to 
electron density fluctuations, as calculated in the next chapter, at a level of 
Sh (J> (Q) =• 1()~ 18 rad 2 /Hz for a carrier density at threshold of 10 1B cm * 3 . Then 
from equation (67), the minimum linewidth for the laser is calculated to be 
120 Hz. This is 30 times less than the lowest reported linewidth of 3.6 kHz 
for a single- mode laser diode. In other semiconductor lasers, the minimum 
observed linewidth is 2 to 4 orders of magnitude larger. This estimate of the 
minimum linewidth implies noise mechanisms other than electron-density 
fluctuations may be dominant in the lasers tested, but that the lowest re- 
ported value to date is not far from the minimum due to electron density 
fluctuations. 

Various mechanisms have been advanced to explain the origin of the 
power-independent linewidth and linewidth rebroadening, including elec- 
tron number-density fluctuations 23], spatial hole burning [{>]. and spatially- 
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dependent temperature and carrier fluctuations 24] These or other mech- 
anisms, such as injection-current fluctuations, may be operative in various 
combinations in any particular laser. Whatever the source, multiplicative 
phase and gain fluctuations in the active medium will cause a minimum 
linewidth for the output field emission spectrum due to the fact that the ac- 
tive medium experiences thermodynamic fluc tuations that alter the refractive 
index [12]. Multiplicative optical phase fluctuations in the active medium 
due to temperature and density fluctuations arise from the same consid- 
erations responsible for Johnson noise in electronic components, Therefore, 
multiplicative noise represents a fundamental limit to the minimum linewidth 
achievable in a laser. When this limit is reached, further increases in gain, 
or decreases in cavity losses, will not afford decreased linewidth. In fact, the 
linewidth is predicted to rebroaden if the net. gain is increased past the point 
that, the large-angle regime is entered. Also, it is noted that the total effect 
of independent multiplicative noise sources is diminished as more sources are 
added in parallel. This means that the effects of multiplicative noise will 
be more pronounced in structures with small transverse dimensions, such as 
semiconductor lasers and optical fibers, as compared to bulk optical systems. 

3.4 Conclusion 

An expression relating the effect of the multiplicative fluctuations to the 
linewidth of a resonant amplifier was derived using the results of the previ- 
ous small-angle analysis. The amount of linewidth enhancement is calculated 
in the small-angle regime from the power spectral density of phase fluefun- 
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tions obtained from the small-angle analysis. The srnall-angle fluctuations 
are found to lead to an output-power-independent Lorentziau power spec- 
trum for the resonant, amplifier. The power spectrum in the large-angle 
regime broadens with increasing net gain as the numerical results of the last 
chapter illustrated. The onset of the large-angle regime is derived from the 
net round-trip gain and the multiplicative noise variance. For the case of 
multiplicative phase noise, it is shown that this effect produces a minimum 
linewidth of the resonant amplifier, and qualitatively explains the phenomena 
of the power-indenpendent linewidth and linewidth rebroadening in single- 
frequency semiconductor lasers. 

4 References 
References 

[1] R. Drever. et a)., Applied Physics B Pliotophysio And Laser Chemistry 
31, 97-105 (1933). 

[2] C H. Henry. IEEE Journal of Quantum Ele ctronics Qli-18, 259-264 
(1982). C. H. Hei.ry. Journal of Lightwave Technology LT-4, 298-311 
(1986). C. H. Henry. IEEE Journal of Quantum Electronics QE-19, 1391- 
1397 (1983). 

[3 K. V'ahala, A. Yariv, IEEE Journal of Quantum Elect ronics QE-19, 1096- 
NO 1 (1983). K. Vahrila, A. Yariv. FEEE Journal of Quantum Elect i onics 
QE 19. 1102-1109 (1983). 

fr r > 7 


(4j A. E. Siegman, Lasers. (University Science, Mill Valley, CA, 1986). 

[5] H. Yauiazaki, M. Yaraaguchi, M. Kitarauia. IEEE Photonic Tech. Lett. 
6. 341-343 (1991). 

[6] H. Olesen. B. Trornborg, H. E. Lassen. X. Pan, Electron. Lett. 28, 444 
(1992). 

[7] D. L. Macfarlane. E. M. Dowling, Journal Of The Optical Society Of 
America A Optics Image Science And Vision 11. 236-245 (1994). 

[8] A. D. Poularikas, S. Seely, Signals and Systems. (PWS Engineering. 
Boston. 1985). 

[9] .1. W. Goodman. Statistical Optics. (Wiley, New York. 1985). 

40] K. Petermann. Laser Diode Modulation and Noise. T. Okoshi, Ed.. Ad- 
vances in Optoelectronics (Kluwcr, Dordrecht, 1991), vol. 3. 

[11] D. Welford. A. Mooradian. Appl. Phys. Lett. 40. 560-562 (1982). 

[12] W. H. Glenn. IEEE Journal of Quantum Electronics 25, 1218-1224 
(1989). 

[13] J. M. Deiitsch, Physical Review E 48. R4179-R4182 (1993). 

[14] M. R. Young. S. Singh, Physical Review A 38. 238-244 (1988). M. R. 
Young. S. Singh. Opt Lett. 13. 21-23 (1988). 


B-S8 


[15] S. Zhu, Physical Review A 41, 1689-1694 (1990), S. Zhu. Physical Review 
A 45, 3210-3215 (1992), S. Zhu. Physical Review A 45. 8148-8153 (1992), 
S. Zhu, Physical Review A 47, 2405-2408 (1993). 

[16] H. B. Callen, R. F. Greene. Physical Review 86, pp. 702-710, (1952). 

[17] M. Be, i, E. Wolf, Principles of Optics. (Pcrgamon Press, Oxford, 1980). 

[18] J. A. Barnes, D. W. Allan, Proc. IEEE 34, 176-178 (1966). 

[19] A. L. Schawlow, C H. Townes, Physical Review 1.12, 1940-1949 (1958). 

[20] E. I. Gordon, Bell System Technical Journal 43, 507-539 (1964). 

[21] M. Sargent, M. O. Scully, W. E. Lamb, Laser Physics. (Addison- Wesley, 
Reading, Mass., 1974). 

[22] K. Kikuchi, IEEE Jour. Quant. Elec. 25, 684-688 (1989). 

[23] D. Welford, A. Mooradian, Appl. Phvs. Lett. 40, 560-562 (1982). 

[ 24] R. J. Lang. K. J. Vahala. A. Yariv. IEEE Jour. Quant.. Elec. QE-21. 
443-451 (1985). 


B-'o 



B - 60 










B- 61 


Ra(E) 



B 


62 


Ra(E 









B - 64 


(3) *3 



Figure - 6. Output Field In Complex Plane for Net Gain, 1^“ 0.999 
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Figure - 18, Relative Probability vs. Intensity, Kc - 0.95 
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Figure - 23. Power Spectrum vs. Frequency, Ko = 0.95 
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Figure - 24. Power Spectrum vs. Frequency, Ko » 0.99 
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OF 

ROME LABORATORY 


Mission. The mission of Rome Laboratory is to advance the science and 
technologies of command, control, communications and intelligence and to 
transition them into systems to meet customer needs. To achieve this, 
Rome Lab: 


a. Conducts vigorous research, development and test programs in all 
applicable technologies; 

b. Transitions technology to current and future systems to improve 
operational capability, readiness, and supportabHity; 

c. Provides a full range of technical support to Air Force Materiel 
Command product centers and other Air Force organizations; 

d. Promotes transfer of technology to the private sector; 

e. Maintains leading edge technological expertise in the areas of 
surveillance, communications, command and control, intelligence, reliability 
science, electro-magnetic technology, photonics, signal processing, and 
computational science. 


The thrust areas of technic* tipetence Include; Surveillance, 
Communications, Command *jid Control, Intelligence, Signal Processing, 
Computer Science and Technology, Electromagnetic Technology, 
Photonics and Reliability Sciences. 


